Coherent-incoherent transition in a Cooper-pair-box coupled to a quantum oscillator: 

an equilibrium approach 
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Temperature effect on quantum tunneling in a Cooper-pair-box coupled to a quantum oscillator is 
studied by both numerical and analytical calculations. It is found that, in strong coupling regions, 
coherent tunneling of a Cooper-pair-box can be destroyed by its coupling to a quantum oscillator 
and the tunneling becomes thermally activated as the temperature rises, leading to failure of the 
Cooper-pair-box. The transition temperature between the coherent tunneling and thermo-activated 
hopping is determined and physics analysis based on small polaron theory is also provided. 

PACS numbers: 85.85.-|-j, 85.25.Cp, 71.38.-k 
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I. INTRODUCTION 

Since the experimental realization of the single- 
electron transistor (SET) in 1987^ the SET has 
become an important element in both extremely 
precise measuremenb^'^'^ and quantum information 
processingiiiSiS The development of the so-called radio- 
frequency SET electrometer is considered as the electro- 
static "dual" of the well known superconducting quantum 
interference devices (SQUIDs) and has shown its power 
of fast and ultrasensitive measurement j On the other 
hand, the direct observation of the coherent oscillation 
in superconducting SET (SSET, or the so-called Cooper- 
pair-box (CPB)) shows the SSET can be served as a phys- 
ical realization of a coherent two-state system)^!^^ and 
since then the SSET or the SSET-based circuit has been 
considered as a candidate of a qubit in a solid state device 
and a read-out device of a qubit. '''^'^ 

The high frequency (up to 1 GHz) and high qual- 
ity factor (Q ~ 10'^) nano-mechanical resonator (NR) 
was firstly fabricated from bulk silicon crystal in 
1996fi^ and therefore it brought a lot of interest 
in demonstrating the quantum nature of this small 
mechanical devicCi ^^'^^i^^'^^i^^ A nano-mechanical res- 
onator capacitively coupled to a SSET forms the so- 
called nano-electromechanical system (NEMS), a sys- 
tem shows promises for fast and ultrasensitive force 
microscopyi^i^ Especially, displacement-detection ap- 
proaching the quantum limit by this system has been 
demonstrated.— From a view point of theoretical study, 
NEMS can be interesting in their own right. In fact, 
NEMS provides a physical realization of a simple quan- 
tum system: a two-level system (TLS) coupled to a single 
phonon mode, which has basic relation to some impor- 
tant models in solid state physics and quantum optics. 
First of all, it represents an oversimplified spin-boson 
model,— 1^ i.e., the spin-boson model in a single- phonon 
mode case, also it can be considered as the Einstein model 
in a two-site problem in the exciton-phonon (or polaron- 
phonon) system.— By taking the analog of the resonator 
as the single mode cavity, it is the Jaynes-Cummings 
model in quantum optics.— iSii^^iS^ 

To be a physical realization of a TLS, the SSET needs 



to work in severe constrains. It is well known that both 
thermal and quantum fluctuations from the environment 
can destroy the Coulomb blockade of tunneling and lead 
to the failure of the SET.^^'^^ Additional conditions are 
needed for a SSETfSI and one important factor is to 
maintain the coherent tunneling. The study on polaron- 
phonon system shown that the coherent motion of the 
electron can be destroyed by its coupling to the phonon 
and a transition from band motion to hopping motion 
happens as the temperature risesi^ By analogy, it is im- 
portant to know if the coherent tunneling of the SSET in 
NEMS can be destroyed by its coupling to a NR, leading 
to failure of the SSET. We shall state such a phonon- 
induced transition as coherent-incoherent transition and 
this is the main interest of the present paper. 

The coherent-incoherent transition in NEMS will be 
studied by both numerical and analytical analyses. It is 
shown that the coherent-incoherent transition exists for 
some frequency ranges with strong coupling parameters. 
The rest of the present paper is organized as follows. In 
Sec. II, we explain the model and the way to trace out 
the environment, the results from adiabatic approxima- 
tion is also discussed. The coherent-incoherent transition 
is studied in Sec. HI and conclusion and discussion are 
presented in Sec. IV. 



II. THE MODEL AND EXPLANATION OF THE 
APPROXIMATION 



The Hamiltonian of the NEMS, i.e. 
given by (setting h = 



a TLS + a NR, is 



CgVg/{2e) with Cg and 



Hf) ^eo^z - ^AoCTa: + ^a)a + \(T^{a) + a), (1) 

where eq = Ec{l ~ '^ng) and Ec is the charging en- 
ergy of a Cooper pair 

V, 

Ao = Ej = Ic/{2e) and /, 
Josephson junction, ai {i - 
ces and (a) is the creation (annihilation) operator of 
the phonon mode with energy O, while A is the coupling 
parameter. The above Hamiltonian indicates that our 



g are the gate capacity and gate voltage, respectively. 



is the critical current of the 
y, z) are the Pauli matri- 
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main interest is the case when the SSET is closed to the 
degeneracy point. It should be noted that the system 
described by the above Hamiltonian is not a thermody- 
namical system, a true thermodynamical system should 
include the environment, which is modeled as a collection 
of harmonic oscillators .-^^^ The whole Hamiltonian can 
be written as 



H — Hq + He 



k 



(2) 



where Hint represents the coupling between the environ- 
ment and NEMS. In the present paper, we shall assume 
that the effect of the environment is just to keep the 
NEMS in equilibrium. Here we apply the concept of re- 
duced density matrix to explain the approximation. The 
density matrix of the whole system is 



Pa ® Pe, 



(3) 



where po and pe are the density matrix of the NEMS and 
environment, respectively. The expectation value of any 
operator Q that acts only on the variables of the NEMS 
can be found by'^- 



(Q) = TTpr{(j)){(j)\Q\(j)), Pr = TrbathP, 



(4) 



where pr is the reduced density matrix and "Trf,at/i" 
means the trace operation over the environment. Our 
approximation means that the above expression can be 
rewritten as 



(Q) 



{4'i\Q\4'i 



Ho\cl>,)^E,\<j>,), (5) 



with Z — J2i 6 is the partition fimction of the 
NEMS and /3 = 1/T (setting fcs = 1). In previous 
treatments,—'^ the effect of the environment is just to 
keep the NR in equilibrium. The present treatment, tak- 
ing a small step further, suggests that the NEMS is in 
equilibrium with the environmental bath. 

By taking this approximation, thermodynamical prop- 
erties of the NEMS can be known, provided that the 
eigenvalue problem of Hq in Eq. (1) can be solved. The 
eigenvalues of Hq can be obtain by the standard numer- 
ical diagonalization technique. In a practical treatment, 
the number of eigenvalues treated must be finite, we 
therefore truncate the eigenspectrum by just taking the 
first smallest N eigenvalues, i.e., 

1 1 ^ 

i i—1 

(6) 

while N is determined by the 



with Z ~ J2^=l 

following condition 



P{En - El) > L, 



(7) 



where L is a fixed number to control the calculation er- 
ror and in the present paper we have i > 20. Accord- 
ingly, one can find the interested physical quantities from 



the numerical results, for example, the free energy of the 
NEMS is 



N 



F(T)^-Tln^, 



(8) 



from which all the thermodynamical quantities of the 
NEMS can be found. In our numerical calculation, we 
have randomly checked the dependence of the results on 
the value of L. We found no observable dependence on L 
for L > 20 by extending L to 50 or even 100. It should be 
noted that, for calculating the temperature dependence 
of tunneling splitting, Eq. (5) should be modified as 



1 ^ 

A(T)c.-^e-^^^|(0,|<7.|0.)|, 



(9) 



this is because the sign of {4>i\(jx\4>i) is not relevant in the 
calculation of A(r). 

Before going to numerical results, we shall first discuss 
the adiabatic approximation in the case of eq = which 
was shown to be a good approximation for Ao/^i < iJ^ 
In the adiabatic approximation, the eigenvalues of Hq are 
characterized by pairs of energies that are given byi^ 



A 



An 



E^{n) ^nn± Aain) - -, A,(n) = 

(10) 

where \n±) = exp{^(A/r2)(a+ — a)}|n) are the displaced- 
oscillator states and n = 0, 1, 2, • ■ • . Taking these eigen- 
values, one can find the temperature dependence of the 
tunneling splitting according to Eq. (9) and the result is 

oo 

A,(r) ^—Y. Aa(n)e~'5"" 2 cosh(/3A„(n)), (11) 

n=0 

where Za = Er?=o e^'^"" 2 cosh(/3Aa(n)). In the limit of 
Aa('T-) <C r^, it can be found that 



Aa(T) ~ ^ Aa{n)pth{n), pth{r 



(12) 

a result that was used in the previous treatments ji^iii 



III. COHERENT-INCOHERENT TRANSITION 

To perform numerical diagonalization, one needs to 
represent the Hamiltonian Hq with a suitable basis. In 
the treatment of Ref. [HI, the basis of | T: i) ® l^-i) is ap- 
plied, where | 1) and | |) are the eigenstatcs of ct^. How- 
ever, with this basis, the elements of the Hamiltonian ma- 
trix are flooded with the overlap terms between different 
states \n±) as shown in the Eq. (8) of Ref.fll. Such over- 
lap terms are rather complicated to compute and there- 
fore the generation of the Hamiltonian matrix will be 
slow. In our treatment, for convenient, we use a very 
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simple basis that consist of the eigenstates of cr^ and a^a 
to represent the Hamiltonian, i.e., a basis of | f, |) ® 
Hence, we represent the TLS terms with | "f , ].) and the 
NR terms with \n) and do a tensor product, the matrix 
of Ho is generated. Of course, both schemes give the 
same results. The diagonalization can give the eigenval- 
ues of Hq, i.e., {Ei}, and the corresponding eigenvectors 
{4>i}, from which the tunneling splitting {(f>i\c7x\(l)i) can 
be found. Noted that, to reach the accuracy condition 
(7), the matrix size is f2-dependent, the lower the fre- 
quency, the larger the matrix size. In the following, we 
shall firstly concentrate our interest on the unbiased case 
of eo = and the biased case will be discussed in the 
end of this section. The eigenvalues in the case of eo = 
we found are in good agreement with the results in Ref. 
[T6l . i.e., the results of the adiabatic approximation given 
in Eq. (10) work well for fl/ Aq > 1. The corresponding 
tunneling splitting of the eigenstates for A/f2 = 0.5 and 
some typical values of Aq/JI are shown in Fig. 1, where 
the results of the adiabatic approximation are also given 
for comparison. 

It can be found that the sketches of both curves show 
good agreement for Ao/51 < 1 and the higher the ex- 
cited states the better the agreement. However, a closed 
look at the curves shows that obvious discrepancy ap- 
pears for both the ground state and the first excited state 
even when Aq/V. = 1/3. Figure 2 shows the details of 
the discrepancy. We found that the discrepancy begins 
to appear even when Aq/H = 0.1 and becomes obvious 
when Ag/n = 1/3. The present results show that, if 
one compares the eigenvalues of the adiabatic approxi- 
mation with the numerical diagonalization, the adiabatic 
approximation seems to work pretty well for Ao/fi < 1; 
however, if one measures the tunneling splitting, the adi- 
abatic approximation works well only for Aq/H, < 0.1. 

After finding out the eigenvalues and the correspond- 
ing tunneling splitting, one can calculate temperature de- 
pendence of the tunneling splitting according to Eq. (8). 
For all the frequencies we studied (down to Ao/57 = 50), 
it is found that the tunneling splitting decreases with 
temperature at the beginning in low temperature regions. 
However, when the coupling strength A/il increases to 
some critical value, A(r) has a "upturn" at some tem- 
perature Tt, which is A/f2-dependent. Typical results are 
shown in Fig. 3. 

The curves shown in Fig. 3 indicate that, for a given 
value of Ao/n and in the low temperature region, the 
main effect of the NR on quantum tunneling is to de- 
creases the tunneling splitting in the weak coupling 
regions. The situation is more or less the same as 
that in the polaron-phonon (or exciton-phonon) system, 
say, coupling to phonons makes the electron become a 
"dressed" one and lowers the hopping rate of the electron. 
However, as the coupling strength increases, the effect of 
the NR on quantum tunneling changes as temperature 
increases to Tt, e.g., the tunneling splitting is enhanced 
by a NR as temperature increases further. The result is 
reminiscent of the transition from Bloch-type band mo- 



tion to phonon-activated hopping motion in the polaron- 
phonon (or exciton-phonon) system. Moreover, Ao(T) 
obtained from Eq. (11) does not show such "upturn" for 
the frequency regions we studied. This implies that the 
"up-turn" is a non-adiabatic effect. 

As we have mentioned in Sec. I, the present model can 
be considered as an analog to the polaron-phonon system 
of Einstein modeL^^ We therefore employ the small po- 
laron theory to analyze the temperature dependence of 
the tunneling splitting. At the beginning of the low tem- 
perature regions, the expectation value of phonon num- 
ber is small, the main contribution to the tunneling split- 
ting is the so-called diagonal transitions. The diagonal 
transition rate can be found by following the way given 
in Ref. 29. Firstly, a canonical transformation is applied 
to the Hamiltonian in Eq.(2), i.e., 

H' = e^He-'^ ^H^ + V + Y,LOkb{h + HLt. (13) 

k 

where S — {\/Vl)<Jz{a^ ~ a), 

Ho = |<T, + r!ata-^, (14) 

V=^^^^e^^/ma-a^) +h.C.), (15) 

^'int ~ s^Hinte~^ and d'± — ± id-y. The diagonal 
transition rate is given by 

w,{T) = Ty^^^{i\V\i)^^e-'-, (16) 

where St = {2X/nf{n + 1/2) and n = (e'''^^ - 1)-^ 
is the expectation value of phonon number. Here, we 
have employed the approximation elucidated in Sec. II, 
i.e., the effect of J^k^kblH + Hint is just to keep the 
NEMS in equilibrium. Wd and hence the tunneling split- 
ting will decrease with increasing temperature. On the 
other hand, the main contribution to tunneling splitting, 
at higher temperature, is the so-called non-diagonal tran- 
sitions. One can calculate the correlation function in the 
way given in Ref. [29l. e.g., 

wit)^Trfy^^mvit)v{om-mn)\% (i?) 

with V{t) = e^^ot^g-^Wot and the result is 

Wit) = (Ao/2)2e-2^-K(*) - /o(e)], (18) 
where /o(a;) is a Bessel function and 

tp{t) = ecos[n{t + iph/2)], e = 2{2\/nf [n{n + 1)]^/^ 

(19) 

It is interesting to note that, in the present case, the 
factor /o(e) still survives after taking thermodynamical 
limit, a situation that is different from the small polaron 
theory of the Einstein model. This is because only one 
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FIG. 1; Tunneling splitting 
for eigenstates of Ho in the 
case of X/Q = 0.5 and eo = 

for some typical values of 
Ao/Q by numerical calculation 
(shown as triangles). The re- 
sults of the adiabatic approxi- 
mation (shown as "*") are also 
shown for comparison. The 
sketches of both curves show 
good agreement for Ao/f2 < 

1 and the higher the excited 
states the better the agree- 
ment. However, obvious dis- 
crepancy appears for both the 
ground state and the first ex- 
cited state even when Ao/fl = 
1/3. The results of the adia- 
batic approximation show large 
discrepancy with the numerical 
results in low frequency region 
with Ao/fl > 1. 
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FIG. 2: Comparison of the tunneling splitting of both the 
ground state and the first excited state by numerical calcu- 
lation with adiabatic approximation in the case of eo = for 
some typical values of Ao/f2. It should be noted that the re- 
sult by adiabatic approximation is Ao/il- independent, which 
is shown as solid and dash lines (dash line is the ground state). 
Numerical results are shown in different patterns for different 
values of Ao/fi (the hollow patterns are for the ground state). 



phonon mode coupled to the TLS in the present case and 
all the other modes just serve as the bath. This result 
can also help to get rid of the delta function problem 
shown in Ref. ^2^. The non-diagonal transition rate can 



be found by using the saddle-point integration 

/oo 
W{t)dt ~ (Ao/2)2(^/7)i/2e-2ST+.^ (20) 
-oo 

where 7 — X^[n{n + 1)]^^^. It can be easily checked that 
Wn shows opposite temperature dependence to Wd and 
the transition temperature Tt is determined by Wd = w„ , 
which leads to 



Ao 



l/2^e{Tt)-ST{Tt) ^ 



1. 



(21) 



The above equation can be solved numerically and the 
transition temperatures obtained are shown in Fig. 4, 
where the transition temperatures determined from A(T) 
by numerical calculation (i.e., curves shown in Fig. 3) are 
also presented for comparison. 

It can be seen from Fig. 4 that the transition tem- 
peratures from numerical calculations are always lower 
than the analytical ones. In high frequency region 
with Aq/D, < 1, both obtained transition tempera- 
tures show similar A/il dependence, i.e., Tt increases as 
A/r2 decreases. However, the transition temperature ob- 
tained from analytical calculation shows opposite cou- 
pling parameter dependence in low frequency region with 
Ao/fi < 1, a result that is in confliction with an intuitive 
picture. We notice that, as shown in the inset of Fig. 4, 
analytical calculation predicts a transition for coupling 
strength as low as X/ft ^ 0.2 for Ao/fi < 1/3; however, 
numerical result shows that there is no transitions for 
coupling strength X/fl lower than 0.7. Figure 5 draws 
the the transition boundary as a function of X/fl and 
Ao/n obtained by the numerical calculation. 
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FIG. 3: Temperature depen- 
dence on the tunneling split- 
ting in the CEise of eo = for 
some typical values of Ao/f2. 
At the beginning, A(T) de- 
creases with temperature, but 
as the coupling parameter A/fi 
increases, A(T) shows a "up- 
turn" at some temperature Tt, 
which is A/f2-dependent. 
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We believe that the discrepancy mainly comes from 
the approximation made in the analytical calculation of 
the diagonal transition rate. It can be shown that 

oo 

WA (T) = ^ A, (n)pth (n) ~ A„ (T) , (22) 

n=0 

which indicates that effectively the diagonal transition 
rate is found by adiabatic approximation, an approxi- 
mation is good for finding the tunneling splitting only 
when Ao/ri < 0.1. Figure 6 shows the details of how the 
discrepancy comes from the adiabatic approximation. 

It turns out that, as the temperature increases, the 
descending rate of Wd{T)/wci{0) is much lower than 
A(T)/A(0) from the numerical calculation. Such a re- 
sult has two consequences. The first one is the transition 
temperature by analytical calculation is higher than the 
numerical result since a slower descending Wd{T)/wd{0) 
will lead to a higher crosspoint with a given w„(T); Still 
another is the extremely high transition temperature in 
the weak coupling regions shown in the inset of Fig. 4. 
It is found that, in the weak coupling regions, the slow 
descending Wd{T)/wd{0) can still survive to some high 
temperature, at where A(T) from the numerical calcula- 
tion has died out, bringing about an artifact high tran- 
sition temperature. On the other hand, As one can see 
from Fig. 1 and Fig. 2, tunneling splitting from the adi- 
abatic approximation shows large discrepancy with the 
numerical result when Aq/CI > 1. This implies that 
Wd{T) is a bad approximation to calculate the diagonal 
transition rate in low frequency regions with Ao/fi > 1. 



Numerical analysis indicates that WdiT) shows different 
(A/r2)-dependence from A(T), leading to different A/f2- 
dependence of Ti. In other words, the break-down of the 
adiabatic approximation in the low frequency region with 
Ao/0 > 1 leads to an abnormal A/fi-dependence of the 
transition temperature in the low frequency regions. 



Now we turn to see the effect of bias on the coherent- 
incoherent transition. In the case of eo 7^ 0, intuitively, 
the hopping between the TLS needs the assistance of the 
phonon unless the bias is very small comparing with Aq, 
i.e., the bias can be overcome by the quantum fluctua- 
tion. Accordingly, the diagonal transition is still possible 
and hence the coherent-incoherent transition is expected 
to survive only when co/Aq <C 1. Nevertheless, the ap- 
pearance of the non-zero bias will weaken the diagonal 
transition rate. As eo increases to some value, at where 
tunneling becomes impossible without the assistance of 
phonon, then the diagonal transition makes no contribu- 
tion to tunneling and coherent-incoherent transition dis- 
appears. Numerical calculation in the case of eo 7^ is 
the same as eo = and some typical results obtained by 
the numerical calculation are shown in Fig. 7. As eg/Ao 
increases, the transition temperature decreases and the 
transition disappears when eo/Ao reaches some critical 
value which is a function of A/f2 and eo/0. Obviously, 
the numerical result is in accord with the analysis pre- 
sented above. 
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FIG. 4: Transition temper- 
atures obtained by analytical 
(from Eq.(21)) and numerical 
calculations in the case of eo = 
for some typical values of 
Ao/n. The analytical results 
are always higher than the nu- 
merical ones. For Ao/f2 < 
1, both transition temperatures 
show similar A/fi dependence, 
i.e., Tt decreases with increas- 
ing A/Jl. However, opposite 
A/n dependence is found for 
A/n > 1, showing that the ana- 
lytical results are invalid in the 
low frequency regions. The in- 
set of the lefttop figure shows 
the whole transition tempera- 
ture curve obtained by the ana- 
lytical calculation. 
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FIG. 5: Transition boundary obtained by the numerical cal- 
culation in the case of eo = 0. The larger the Ao/fi (i.e., 
the lower frequency for a given Ao), the larger the coupling 
strength A/fi is needed for the transition to happens, indi- 
cating the SSET is more stable when coupling to a lower fre- 
quency NR. 




FIG. 6: Illustration of how the discrepancy comes from 
Wd{T) obtained by the adiabatic approximation in the case 
of €o = 0, Ao/Q = 1/3, and A/fi = 1. As the temperature 
increases, the descending rate of Wd{T) /wdiO) is much lower 
than A(r)/A(0) from the numerical calculation, leading to 
a higher transition temperature. The "up-turn" of Wn{T) in 
the low temperature regions is due to the factor 7~^^^ which 
diverges as T — » 0. 



IV. CONCLUSION AND DISCUSSION 

In conclusion, we have presented study on the effect 

of a nanomcchanical resonator on quantum tunneling in 
a Cooper-pair-box at T 7^ 0. We found that the co- 



herent tunnehng of a Cooper-pair-box can be destroyed 
by the coupling to a NR at a temperature much lower 
than the coupling energy Ej of the Josephson junction. 
The present analysis shows that, for a NEMS to work 
well, one additional condition, e.g., T Tt, is needed. 
The coherent-incoherent transition boundary as a func- 
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FIG. 7: The efltect of bias on 
the coherent-incoherent transi- 
tion for some typical values of 
X/Q, and Ao/Q,. The transition 
temperature decreases with in- 
creasing eo /O. and the transition 
disappears when eo/Sl reaches 
some critical value which de- 
pends on both and Ao/Q. 



tion of Ao/n and A/fi is calculated. It turns out that 
the transition happens only for X/fl > 0.7, and the 
lower the Ao/fi, the larger the A/il is needed for the 
transition. We also found that the transition tempera- 
ture is a monotonic descending function of X/fl. The 
present result shows that experimental observation on the 
coherent-incoherent transition is still impossible since the 
corresponding coupling parameter cannot be achieved in 
the present stage. Taking some typical experimental 
parameterssi^ Aq = AfieV and hil = 1.2/j,eV, the cou- 
pling parameter needed for the transition is A/f2 ~ 1.4, 
which is larger than the strong coupling limit possible to 
achieve in experiment [X/hvo ~ To see the transi- 
tion in the coupling parameters can be achieved in exper- 
iment, the result shown in Fig. 5 tells Ao/Sl should lower 
than 1, i.e., the frequency of the NR should be as high 
as IGHz, which is almost the limitation in the present 
stage. Nevertheless, it is believed that, in a real system, 
the transition can be seen in a lower coupling parame- 
ter regions since the coupling of the NR to the environ- 
ment can help the transition to happen. It is also found 
that other thermodynamical quantities, like specific heat, 
vary smoothly over the transition point, showing that the 
transition is not a thermodynamical transition. 

Coherent-incoherent transition is an important issue in 
the spin-boson model and analysis on the transition at 
T ^ was also provided in Ref. 2lL- However, the start- 
ing point in Ref. [2ll is different from the present analysis. 
In the present analysis, the key element is the tempera- 



ture dependence on the tunneling splitting, while in the 
previous analysis, it is the relaxation behavior, i.e., the 
time dependence on transition rate P{t).— As a mat- 
ter of fact, the present analytic calculation is similar to 
the calculation provided in Sec. Ill D of Ref. 21, and 
accordingly the transition rate we found here is corre- 
sponding to r (or 1/(2t)), which was predicted to have a 
monotonic temperature dependence in their analysis. It 
seems that this result is suitable for the case with large 
bias while the "up-turn" of A(r) for small (or zero) bias 
cannot be explained. Nevertheless, the present result 
can help to understand the coherent-incoherent transi- 
tion at r 7^ in the spin-boson model with small bias. 
To the first order approximation (i.e., omitting the co- 
operative effect between different phonon modes), the 
combined contribution of all the phonon modes with a 
frequency-dependent weight for < lu < lOc is approxi- 
mately the contribution of the whole bath, accordingly a 
coherent-incoherent transition is expected to exist in the 
spin-boson model since the coupling of the TLS to all 
the phonon modes can lead to the transition when the 
coupling strength exceeds some value. 
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